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Time-reversible symmetry arises frequently in many studies in classical and quantum mechanics. A 



C/3 

q 

.12 Abstract 

We present an algorithm for finding all time-reversible systems within a given family of 2-dim 
systems of ODE's whose right-hand sides are polynomials. We also study an interconnection of 
time-reversibility and invariants of a subgroup of SX(2,C). 

> 

oo ! 1 Introduction 

general introductory account and a survey of state of the art with regards to time-reversible symmetry 
C ■ in dynamical systems theory and its applications to some physical problems can be found in [8]. 

Another recent stream of research is devoted to investigation of isochronicity and linearizability of 
| time-reversible polynomial systems (see, e.g., [TJ [TUJ, [TT] and references therein). 

In the context of dynamical systems described by ODE's of the form 

§ = *■(■) (»€«), (1) 

where F : Q i— > Tf2 is a vector field and O is a manifold, by time-reversible symmetry we mean an 
invertible map R : £1 >—> J7, such that 

^ - < 2 > 

In this paper we investigate probably the simplest time-reversible symmetries, namely, we limit 
our study to the case of 2-dim systems, that is, z = (x, y) G C 2 (or M 2 ), and to the symmetries of the 
form 

R i x i — > 7y, y i— > 7 _1 x (3) 

with 7 £ C (or M) (note, that R defined by ([3]) is an involution). Using methods of computational 
algebra we show an interconnection of the physically important phenomena of time-reversibility, involu- 
tion and group invariants in dynamical systems described by autonomous polynomial two-dimensional 
systems of ODE's. 
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2 Preliminaries 



We say that a straight line L is an axis of symmetry of a real autonomous two-dimensional system of 
ODE's if as point-sets (ignoring the sense of the parametrization by time t) the orbits of the system 
are symmetric with respect to the line L. There are two types of symmetry of a real system with 
respect to a line L: mirror symmetry, meaning that when the phase portrait is reflected in the line L 
it is unchanged; and time-reversible symmetry, meaning that when the phase portrait is reflected in 
the line L and then the sense of every trajectory is reversed (corresponding to a reversal of time), the 
original phase portrait is obtained. The symmetry ([3]) can be considered as a generalization of the 
second type of reflection to the case of complex two-dimensional systems in the following sense. 
Consider a parametric family of real systems 

u = U(u,v), v = V{u,v). (4) 

Introducing the complex structure on the plane (u, v) by setting x = u + iv we obtain from Q the 
equation 

x = P(x,x), (P = U + iV). (5) 
We add to this equation its complex conjugate to obtain the system 



x = P{x,x), x = P{x,x). 

Let us now consider x as a new variable y and allow the parameters of the second equation to be 
arbitrary. Then ([5]) yields the complex system 

x = P(x,y), y = Q(x,y). (6) 

Let o denote the vector of coefficients of the polynomial P(x, x) in ([5]), arising from the real system 
([3]) by setting x = u + iv. It is easy to see that if a = ±a (meaning that either all the coefficients are 
real or all are pure imaginary), then the li-axis is an axis of symmetry of the real system (j4|) and of 
the corresponding complex differential equation ([5]). Thus, u-axis is an axis of symmetry for ([5j) if 



P{x,x) = -P{x,x) (7) 

(the case a = —a), or if 

P(x,x) = P(x,x) (8) 

(the case a = a). We now observe that if condition (JT]) is satisfied then under the change 
equation ([5]) is transformed to its negative, 

x = -P(x,x), (9) 

and if condition ([8]) holds then ([5]) is unchanged. Thus condition ([8]) means that the system is reversible 
with respect to reflection across the it-axis (i.e., the transformation does not change the system) while 
condition ([7]) corresponds to time-reversibility with respect to the same transformation. 

If the line of reflection is not the ii-axis but a distinct line L then we can apply the rotation 
x\ = e~ lLp x through an appropriate angle ip to make L the u-axis. In the new coordinates we have 

Xl = e - itp P(e iip x 1 ,e- iip x 1 ). 
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By the discussion in the paragraph following this system is time-reversible with respect to the line 
Imxi = if ([7|) holds, meaning that 



e i vp( e ^xi,e-^x 1 ) = -e~^P(e^i, e~ i<p xx). 

Hence, reverting to the variable x, ([5]) is time-reversible when there exists a cp such that 

e 2i ^P{x,x) = -P{e 2i ^x, e~ 2i ^x) . (10) 

We see that in the particular case when in ([2]) F = (P, Q), 7 = e 2i(p , y = x, and Q = P the equality 
([2]) is equivalent to (|10p . Thus, the second kind of reflection with respect to a line described at the 
beginning of the section can be considered as a particular case of symmetry ([3]). 

In fact we have shown that if system © is the complexification of a real system @ and it admits 
a symmetry ([3|) with 70 = e 2l ^° then the line v = u tan(/>o is the line of symmetry of the trajectories 
(as point-sets) of the real system (see [13] for more details). 

3 The set of time-reversible systems 

Direct calculation shows that the system ([6]) is time-reversible with respect to a transformation ([3]) if 
and only if for some 7 

lQ(iy,x/j) = -P(x,y), jQ(x,y) = -P{jy,x/~f) . (11) 

We will limit our study to the case when P(x,y), Q{x,y) in ([6]) are polynomials. Without loss of 
generality we can write such systems in the form 

^ = - E a pq xv +1 y<i = P(x,y), §= E V*V +1 = Q(x, y), (12) 
(p, q )eS ( P ,q)eS 

where S is the set S = {(pj, qj) \pj + qj > 0, j = 1, . . . ,£} C ({—1} U No) x No, and No denotes the set 
of nonnegative integers. The notation fjl2[> simply emphasizes that we take into account only nonzero 
coefficients of the polynomials and the symmetric form of the system (in the sense that if there is 
a term a pq x p+l y q in the first equation of (|12p . then there is also the term b qp x q y p+1 in the second 
equation). Obviously, any polynomial system ([6]) can be embedded into a system of the form (|12|) . 

We will assume that the parameters CL Pjqj , b q . p . (j = 1, ... ,£) are from an infinite field k. Denote 
by (a, b) = {a piqi , . . . , a Piqi , b qePe . . . , b qipi ) the ordered vector of coefficients of system (fT2j) . by E(a, b) 
the parameter space of (|12p (e.g. E(a,b) is C 2i or M 2 ^), and by fc[a,6] the polynomial ring in the 
variables a pq , b qp over the field k. It is clear that there is a one-to-one correspondence between points 
of E(a,b) and systems of the form (|12p . The condition (jlip immediately yields that system (|12p is 
time-reversible if and only if 

b qp = l p ~ q a pq , a pq = b qp ^ v . (13) 

We rewrite f)13|) in the form 

a Pkqk = t k , b qkPk = 7 Pfe " 9fc tfc (14) 

for k = 1, . . . ,£. From a geometrical point of view equations (I14p define a surface in the affine space 
k 3e+1 = (a Piqi , . . . , a Peqe , b qePe , . . . , b qiPl , t\, . . . , tg, 7). Thus the set of all time-reversible systems is the 
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projection of this surface onto k = E(a, b). Nowadays methods to compute such projections are well 
known (see e.g. [5] Chapter 3]). Namely, let 

H = (1 - wji ■ • -7£, a Pkqk - t k , %b qkPk - %tk I k = 1, . . . ,£), (15) 

where j k = j qk ~ p «, % = 1 if Pk - Qk < 0, j k = 1, % = j Pk ~ qk if p k - q k > and (hi,..., h m ) denotes 
the ideal generated by polynomials hi, ... , h m . Then as a direct corollary of Theorem 2 of [5, Chapter 
3, §3] we immediately find the minimal variety which contains the set of all time-reversible systems in 
the family (|12p (we remind that the variety of ideal /, denoted by V(7), is the set of common zeros 
of all polynomials of /; the Zariski closure W of a set W in the affine space k n is the minimal variety 
which contains W). Namely, from the theorem 

K = V(1) where T= k[a,b] C\H, (16) 

that is, the Zariski closure of the set 1Z of all time-reversible systems is the variety of the ideal I. 
Moreover, according to Theorem 2 of [5j Chapter 3, §1] in order to find a generating set for the 
ideal X it is sufficient to compute a Groebner basis for H with respect to an elimination order with 
{w,~f,tk} > {a Pk q k , bg kPk } and take from the output list those polynomials, which depend only on 
a Pfe9fe' bq kPk (k = !,...,€). 

We now will obtain another description of the ideal I, which will link the time-reversibility to 
invariants of a subgroup of SL(2, C). Denote by M. the set of all solutions v = (ui, u^, ■ ■ ■ , 1^21) with 
non-negative components of the equation 

Cl^l + Q2V2 H + Q+m+i H H (2£V2£ = 0, (17) 

where Q = pj - qj for j = !,...,£, Q = qu-j+i - P2£-j+i for j = I + 1, . . . , 21, that is, 

C = (Pi - 9i,P2 - Q2, ■ ■ ■ ,pe - qe,qe -pt,...,qi - pi) 

(we remind that (pj, qj) are from the set S defining system (|12H ) . Obviously, M. is an Abelian monoid 
(semigroup). For v = {yi, . . . , V21) G we denote by [v] the monomial 

a Pi9i a P292 ' ' ' a piqt^Qm^Qt-lPt-l ' ' ' KlPi (^^) 

and by v the involution of the vector v, v = (v 2 e, V2i-\i ■ ■ ■ > We call the ideal 

I s = ([ u ] -[P]\v G M) C k[a,b) 

the Sibirsky ideal of system (fT2l) (Sibirsky [13], [H] studied this ideal for the case of real systems 
embedded in the family (j 12|) ; some studies in this direction were performed also in [H [6] [9]). In the 
case that (fT2l) is time-reversible, using (fT3|) and (|T7|l we see that for ^ G 

[i>]= 7 C->] = M > (19) 

where £ • ^ is the scalar product of £ and i/, that is the left-hand side of (|17p . By (|19p every time- 
reversible system (a, 6) G ^(a, b) belongs to *V(Is). However it is easily seen that the converse is false. 
The following statement proven in [12] gives the exact characterization of the set of time-reversible 
systems. 

1 Groebner basis is a cornerstone of many algorithms of computational algebra (see e.g. [5] for the definitions). 
Mathematica, Maple etc. have routines to compute it. 
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Theorem 1 Let 1Z C E{a, b) be the set of all time-reversible systems in the family then 

(a) K c V(/ s ); 

(b) V(Js-) \TZ= {(a, 6) | 3(p, q) € S such that a pq b qp = but a pg + 6 gp ^ 0}. 

In other words, (b) means that if in a time-reversible system (|12|) a p(? 7^ then b qp 7^ as well. It 
follows from (b) of the statement that the inclusion in (a) is strict, that is 1Z ^ V(/s). 
From (|16p and (a) we have that V(T) C V(/s). We will show that in fact 

Is =1 and both ideals are prime. (20) 



When (|20p is proven, from (|16p and (|20p we obtain the following description of time-reversible 
systems. 

Theorem 2 T/ie variety of the Sibirsky ideal Is is the Zariski closure of the set TZ of all time-reversible 
systems in the family 

We now prove (|20p. Suppose we are given the system 

_ fi(h, ...,t m ) _ f n (ti, ...,t m ) 

•'1— /, 1 \ ' ' ' ' 5 TL /, J. \ ' V I 

9l{tl, ■ ■ ■ , tm) 9n{tl, ■ ■ ■ ,t m ) 

where fj,gj € k[t\, . . . , t m ] for j = 1, . . . , n. Let k(t\, . . . , t m ) denote the ring of rational functions in 
m variable with coefficients in k, and consider the ring homomorphism 



4> : k[xi, ...,x n ,ti,.. -,t m ,w] -> k(h, ...,t 



defined by U -> tj, -> . . .,t m )/gj(ti, . . .,t m ),w -> • • • ,i m ), i = 1, • • • , m, j = 1, . . . , n 

and g = gig 2 ■ ■ ■ g n - Let 

H = (1 - wg,xigi(ti, . . . ,t m ) - f\(ti, . . . ,t m ), . . . ,x n g n (ti, . . . ,t m ) - f n (t\, . . . ,t m )). 

It is not difficult to check that 

H = ker(4>). (22) 

Since k[x±, . . . , x n , t\, . . . , t m , w] is a domain (|22|) yields that H is a prime ideal. 
By (|22p the ideal H defined by (|15p is the kernel of the ring homomorphism 

i/j : k[a,b,h, . . .,t e ,-/,w} — ► fe(7,*i,. ..,t e ) 

defined by a Pkqk 1— > tfc, &g fcPfc 1— ► J Pk ~ qk tk, w \— > 1/(71 ••• 7^) for = !,...,£. We obtain a reduced 
Groebner basis G of A; [a, 6] n by computing a reduced Groebner basis of H using an elimination 
ordering with {a Pjqj , b gjp .} < {to, 7, t 3 -} for all j = 1, . . . ,£, and then intersecting it with k[a, b]. Since 
H is binomial, any reduced Groebner basis G of H also consists of binomials. This shows that H is a 
binomial ideal. 

It is easily seen that any binomial of the form [a] — [a] (a = (ai, . . . , a.21) and [a] is defined by 
(|18p ) is in H (which is equal to ker(^)) yielding [a] — [a] G k[a, b] D H = X and, therefore, Is C 1. We 
prove now that 

X C I S . (23) 
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X is a binomial ideal, thus to verify (|23[) it is sufficient to show that all binomials of X belong also to 
Is- Since X is prime it is sufficient to consider only irreducible binomials. Indeed if [a] — [f3] G X is 
not irreducible, then [a] — [f3] = ([ati] — where [ai] — [fii] is irreducible binomial. It is clear 

from the definition of ip that H contains no monomials, thus, since H is prime, [ai] — [fdi] G H. Let 
now [a] — [f3] be an irreducible binomial in /, that is, supp(a) D supp(/3) = 0. Then 

^([a] - [/?]) = t" 1 • • • t^t^ +1 7 (^-*)^+i • • • t ^ 7 (pi-9i)^+i _ 

ff 31 . . . jfii+Pt+i ~,(pe-qe)0e+i . . . ffoe (pi-qi)/3e+i _. 
l l l l l £ I l l 7 — 

A+fai+fo+foi-i f l 3 e+l 3 i+i^(e+iPe+i+i:i+2Pe+2-+P2iP2i 
L l b 2 ■ ■ ■ l l I 

Thus, V ([a] -[/?]) = if and only if 

ai + a 2 i- i+ i = & + fht-i+l for all i = 1, 2, . . . ,£ (24) 

and 

Ce+i a e+i + Gt+2&e+2 1- (210-21 = Ce+ify+i + (1+2^1+2 • ■ ■ + (21(^21- (25) 

Since supp(a) n supp(/3) = 0, ([2"3| yields that 

A = a M _ i+ i fori = l,...,21 (26) 

Hence /9 = (/?i, . . . , = (a^, . . . , ai) = a. Noting that d = —Q21+1-1 and using (|2"6|) we write (f25|) 
as + C2CK2 • • • + 0°^ + (t+iot+i + 0+2^+2 • ■ ■ + C2^«2^ = 0. That means a satisfies the equation 
(fTT|) and, hence, a G AI. This completes the proof of (f20|) . 

Theorem 2 of [H Chapter 3, §1] provides the following algorithm for computing a generating set 
for the ideal X and, therefore, for the ideal Is- 

• Compute a Groebner basis Gh for H defined by (|15p with respect to any elimination order with 
{w,j,t k } > {a Pkqk ,b qkPk I & = 1,... ,£}; 

• the set Gh Clk[a,b] is a generating set for I; 

• H = {f, v I [u] — [u] G G} U {ei + ej | j = 2^ — i+ 1; i = 1, . . . , £} is a Hilbert basis for the monoid 
M (here ej is the element (0, . . . , 0, 1, 0, . . . , 0), with the non-zero entry in the i-th position). 

The correctness of the last step of the algorithm follows from Theorem 4 of [7j . 
As an example consider the real system 

u = —v + a\u 2 + CL2UV + 03V 2 , v = u + b\u 2 + b2uv + 63V 2 . (27) 

Using the complexification procedure described in Section 2 we obtain from (|27p the system 

x = ix + a w x 2 + a ixy + a_i 2 y 2 = ix + P 2 (x, y), ^ 
y = -iy + b 10 xy + b 01 y 2 + b 2 -ix 2 = -iy + Q 2 (x,y). 

In the particular case bk s = ak s , y = x, a\ = ((oj + b 2 - a 3 ) + i(b\ — a 2 — &3))/4, a i = (oj + a 3 + 
i{b\ + bs))/2, a_i2 = (ai — 03 — 62 + + «2 — bs))/4 system (|28|) is equivalent to ([2"?]h 
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Computing a Groebner basis of the ideal 



J = (l - W7 4 , aio - *i, &oi - 7*1) «oi - *2, 7^io - h, a-12 ~ h, 7 3fo 2 -1 - £3) 

with respect to the lexicographic order with w > 7 > t\ > t 2 > t% > aio > a oi > a -i2 > ^10 > &01 > 
62,-1 we obtain a list of polynomials. (The list is too long to be presented here, but is easily computed 
using any computer algebra system.) According to step 2 of the algorithm in order to obtain a basis 
of Is we just have to pick up the polynomials that do not depend on w, 7,ii,i2>*3- In fact, there are 
five such polynomials in the list: /1 = a^ -1 -a_i 2 6? , h = ai a 01 - 601610, h = a? a_i 2 - 62,-16^, 
U = aioa-i26? - 0(^62 _i6 i, fe = a? a -i26io - 00162,-1601- Thus, for system (j2EJ) 

Is = {fl, ■ ■ ■ , fs}- 

An algorithm for computing generators of the ideal Is has been also obtained in [7j. To compare 
the efficiency of the algorithms we carried out computations for few families of system (|12p on a PC 
with a 1.7 GHz processor and 1.5 Gb RAM using Mathematica 4. It takes 0.2 seconds CPU time to 
find Is using the algorithm described above and 1.7 s CPU using the algorithm of [7j. Replacing in 
(|28p the polynomials P2 and Q2 by the homogeneous polynomials of the third and fourth degrees we 
find that it takes, respectively, 0.04 and 4 s CPU with the algorithm of the present paper, and 480 
and 33450 s CPU with the algorithm of 0. 

4 Time-reversibility and invariants 

We now show an interconnection of time-reversibility and invariants of a group of transformations of 
the phase space of system (fT2|) . Let G be a matrix group acting on x = (x%, . . . , x n ) and let k be any 
field. We recall that a polynomial /(x) E fc[x] is invariant under G if /(x) = /(j4 • x) for every A £ G. 
Consider the transformations of the phase space of (|12[) 

x' = t]x, y' = r]^ l y (x,y,r)EC, 77 ^ 0). (29) 

The transformations (129j) form a subgroup of SX(2,C). In (x',y r ) coordinates system (112p has the 
form 

= Yl a ( 7 ?)(p,g) x ' p+1 2/ /9 ' v = b( (n)(q,p) x ' q y' p+1 

(p.9)gs (p,<z)eS 
and the coefficients of the transformed system are 

a(v) Pkqk = a Pkqk V Qk - pk , b( V ) gkPk = b qkPk rf>*-*\ (30) 

where k = 1, ...,£. Let U„ denote the transformation (|30p . f/„ is a linear representation of group 
([29]) in C 2 ^. We will write ([30]) in the short form (a(rj) , 6(r/)) = U v (a,b). It is straightforward to see 
that a polynomial /(a, 6) G C[a, 6] is an invariant of the group U„ if and only if each of its terms is an 
invariant. 

The action of on the coefficients ay , bji of the system of differential equations (I12p yields the 
following transformation of the monomial [u\ defined by (|18|) : 

UM = <v);\ qi ■ ■ • a(vT p * m K'nfqm ■ ■ ■ b «U = (31) 

'I u P iqi u Piq£ u 1lVl u q\V\ ~ ' L J" 
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Thus we see that the monomial [u] is invariant under the action of U v if and only if £ • v = 0, i.e., if 
and only if v G M.. 

For h G C 2e let [z/]|/j denote the value of the monomial [v] = o? v \ a V2 ■ ■ ■ b U21 at the vector h, 
that is, [f]\h = h^h 1 ^ ■ ■ ■ h^f ■ Denote by (a, b) the involution of (a, 6), 

(a, 6) = (b QlPl , . . . , , a P; ,j ; , . . . , a Piqi ). (32) 

Let O be an orbit of the group and assume that there is (a* ,6*) G such that (a*, 6*) G O. By 
(|13p and f|30|) the system (a*, 6*) is time-reversible. We prove that then for any system (a , 6°) G O 

the system (a , 6°) is also in O. Let (a , 6°) = U m {a* ,&*). Any invariant [i/] of the group is 
constant on any orbit of the group, (a*, b*) and (a*,b*) belong to the same orbit O of the group 
U v , hence for all invariants [v] of U v M|( .,&*) = MI ( jT£ry Therefore, H|( a *,b*) = [i>]|( a * ,6*) yielding 
M!c/, ;o (a*,fe*) = K la^ (a*, 6*) and, hence, 

[v}\u no (a*,b*) = Mu Vo{a *,b*y ( 33 ) 

where [/^(a* ,6*) is the involution of the vector [/^(a*, 6*). From (|33|) and Theorem 3.2 of [T2] we 
conclude that L^ (a*,6*) belongs to the orbit O and, therefore, the system U m {a* ,&*) = (a , 6°) is 
time-reversible. 

We say that the orbit O of the group U v is invariant under the involution ([32]) if for any (a, 6) G 
the system (6, a) also belongs to O. By our reasoning above, we have proven 

Theorem 3 (a) The set of the orbits ofU^ is divided into two not intersecting subsets: one consists of 
all time-reversible systems and only time-reversible systems, and there are no time-reversible systems 
in the other subset. 

(b) The variety V(/s) is the Zariski closure of all orbits of the group invariant under the 
involution A32\) . 

To finish we note that the method described in the present paper can be applied to compute the set 
of systems time-reversible with respect to some other symmetries, also in higher dimensions. However, 
we believe that the beautiful algebraic structure exhibiting in the parameter space by the symmetry 
([3]) is due to the fact that ([3]) is an involution. Note also that an interesting future project would be 
a generalization of the presented results to the case of higher dimensional systems. 
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